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Abstract 

We investigate the spreading of thin, circular liquid drops of power- 
law rheology. We derive the equation of motion using the thin film ap- 
proximation, construct source-type similarity solutions and compute the 
spreading rate, aparent contact angles and height profiles. In contrast 
with the spreading of newtonian liquids, the contact line paradox does 
not arise for shear thinning fluids. 

In this work wc study the spreading of circular drops of power-law rheology 
fluids, also known as Ostwald-de Waele fluids Pj. The power-law rheology is 
one of the simplest generalizations of the Newtonian one, in which the effective 
viscosity is assumed to be a power law of the local rate of deformation 7 as 
H = m\j\^^^^^ . The values of m and A depend on the physical properties of the 
liquid. When A > 1, the fluid is called shear-thinning and the viscosity tends to 
zero at high strain rates pp. 

The problem of drop spreading has been intensively studied in the last 
decades (see [21 El EI El E])- The motivation is that this class of flows plays 
a very important role in processes such as coating and painting. One reason for 
our study is that while the past work considers only newtonian fluids, most of 
the fluids of technological interest are non-newtonian. 

There is also a theoretical motivation: one modelling difficulty for newtonian 
fluids is the paradox of the contact line, that states that the dissipation of energy 
is unbounded near the advancing contact line of a newtonian fluid. This paradox 
may be solved, for example, by introducing the effect of the intermolecular 
forces, or abandoning the no-slip boundary condition at the substrate. Here 



we show the remarkable fact that for power-law fluids with A > 1 the paradox 
does not arise. Moreover, there are solutions with moving interfaces, and their 
energy dissipation rate is bounded. The qualitative reason lies in the fact that 
the flow creates a layer adjacent to the substrate where the strain rate diverges, 
and thus, for A > 1 the effective viscosity rj tends to zero. In other words, the 
drop self-lubricates. 

Spreadings of power law fluids were first studied experimentally in and 
the one dimensional flow was analyzed in [S] |^. Here we construct for the 
first time radially symmetric spreading solutions, which can be used directly to 
model droplet spreading. 

In power-law fluids, the deviatoric stress tensor is related to the strain 
tensor according to the following constitutive relation [TllTn): 



T,. 



^m\j\'/'-%j, (1) 



where I7I = y 5 J2i j Ty'Tu ^^'^ Ty strain tensor, given by 

dvi dv4 

Here w,; is the fluid velocity field. When A = 1, one recovers a Newtonian fluid. 
If A > 1, the fluid is said to be shear-thinning, which is the most common case 
that includes many polymer solutions. The values of A arc typically between 
1.7 and 6.7 [J- 

In order to derive the equation of motion we use the lubrication approxima- 
tion, which is widely applied in the case of Newtonian flows (for which A = 1) 
(see for instance [HJ |2] and references therein) . We assume that the film is much 
thinner than its horizontal dimension, that the motion is nearly horizontal and 
that the inertial effects are negligible so that the flow is governed by a balance 
between capillary and viscous forces. We use a cylindrical coordinate system 
(r, 6, z) and assume circular symmetry In this system, the velocity of the fluid 
has components {u{r, z, i), 0, v{r, z, t)) and the substrate is at z = 0. We neglect 
systematically the z component of the velocity when it is compared with the r 
component, \u\ » \v\. We also assume that the stresses are mainly due to high 
gradients of the horizontal velocity u in the z direction. Consistently, we assume 
that the components jrz and Trz are much larger than all its other respective 
components. Then, the r-component of the momentum equation is reduced to 

where p is the pressure. Let z — h{r,t) be the fluid free surface. The conserva- 
tion of mass can be written within the lubrication approximation as 

ht + -{rUh)r = 0, (4) 
r 

where U is the horizontal velocity averaged in the z coordinate. We assume that 
the drop size is much smaller than the capillary length. Then we can neglect 
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gravity and the pressure under the free surface may be approximated by 

p{r,t) = 7K + P0 (5) 

where 7 is the surface tension, n is the mean curvature of the free surface and po 
the atmospheric pressure. Under the lubrication approximation the curvature 
can be approximated by 

K ~ —hrr -. (6) 

r 

One can integrate Eq. El with respect to z to obtain the shear stress 

Trz=Priz-h), (7) 

that satisfies the zero-tangential stress condition at the free surface z — h{r,t). 
By using the constitutive relation given by Eq. ^ 

Tr, =m\u,\^/^'^u„ (8) 

and by substituting Trz from Eq. \7\ we can compute the z derivative of the 
horizontal velocity, 

u, = -m-^\pr\^-^Prih-z)^. (9) 

After integration and using the no-slip condition 
u{z = 0) = one obtains 



.^^-H-V^^^^-^j. (10) 

Now, in order to apply the equation for the conservation of the mass Eq. 0] we 
need the averaged velocity U, 

U=j^J^ udz = -ni-^\pr\^-^prj^. (11) 
Finally, using Eq. ^ we obtain the equation of motion 

1 1 r7A\ ,A+2| ,A-1 



ht--^—:[-) irh^+'\KX-'Kr)r=^0, (12) 
r \+2 \m/ 

which forms a system of partial differential equations with Eq|Sl 

Now we solve the problem in which a finite volume V of fluid is initially 
concentrated on a point in a planar surface. As there are no external length- 
scales in this problem, we shall look for similarity solutions of the first kind 

m 

where A and B are constants and /? is the similarity exponent. The argument 
of H is the similarity variable rj — r/At^^ . If the function H{ri) has a zero at 77/ 
then the front of the spreading will be given by 

rf{t)=Aiift''. (14) 
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Eq. El guarantees that the fluid volume V is conserved, i.e. that 

/•'■/(*) 

V ^ 2TTrh{r, t)dr = A^BI = const. 

Jo 

where we define the shape factor, 

fVf 

/ = / 2nr]H{r])dT]. 

JO 

It is convenient to introduce the dimentionless curvature 

K = -H" - — (15) 
V 

such that, from Eq. K(r,t) = K{'q)V/iA'^It'^'^). By substituting Eq. [HI in 
Eq. 1121 we find that the similarity exponent is given by 



7A + 3' 



By writing A'^^^'^ = {^/m)^{V/I)^^^^ and integrating once, we obtain the 
ODE 

-7]H = H^+^\K'\^-^K' (16) 

which has to be integrated with Eq. E| It is clear that for ?7 > and H > 0, 
K' < and then the equation can be written as 

r]H = H^+^{-K'f. (17) 

The initial conditions are 

i?(0) = 1, H'{Q) = 0, K{Q) = kq. (18) 

The first condition represents the scaled height at the center of the drop, the 
second is the radial symmetry condition, and the third is the central value of 
the curvature, which is a free parameter kq that has to be calculated. 

Using Eq. I14l we obtain a simple formula to compute the radius of the drop 
as a function of the volume V, 7 and m: 

rf{t)^S,{V''+\^/m)Hf^'''''\ (19) 

where Sx = ((7A + 3)/(A + 2))1/(7A+3),^^/j(2A+i)/(7A+3) depends only on A. 

The aparent contact angle of the drop is defined as tan0(t) — —hr{ri,t), 
where is the point of maximum slope, hrr{ri,t) — 0. By simple substitution, 

tan 0(t) = Qa1^V(7A+3) (^/^)3A/(7A+3)^-3/(7A+3) (20) 

where Qx = -H'{f]^){{X + 2)/(7A + 3))3/(7a+3) j-a/(7a+3) ^^^^ ^^^^^ 
where H"{rj,) = 0. 
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A 




Vf 


I 


Sx 


Qx 


1.1 


7.64300 


0.73220 


0.83124 


0.86880 


1.814294 


1.5 


4.08651 


1.03906 


1.61381 


0.99651 


1.137746 


2.0 


3.29562 


1.19915 


2.08030 


1.05262 


0.951946 


2.5 


3.00428 


1.28882 


2.34984 


1.08074 


0.878580 


3.0 


2.85609 


1.34757 


2.52782 


1.09766 


0.840258 


4.0 


2.71025 


1.42076 


2.74976 


1.11683 


0.802190 


5.0 


2.64043 


1.46484 


2.88307 


1.12721 


0.784121 



Table 1: Constants defining the selfsimilar solutions. 



We solve Eqs. 1151 and [T7I sub i ect to Eq. 1181 numerically with a fourth order 
Runge-Kutta scheme, with a variable stepsize equal to iJ/10000. In Fig. 1 we 
show several height profiles for different values of kq with A = 5/2 fixed. For 
kq < 3.00428 the solution is always positive: this solution does not represent a 
drop. For kq > 3.00428, the drop profile tends to zero linearly at 77 = rjf. For a 
critical value of kq « 3.00428 the profile tends to zero with zero contact angle. 
This case represents a drop of a liquid that wets the surface. A similar picture 
is obtained for all A > 1. In general, there is a solution with zero contact angle 
for a critical value of kq = ka and a continuum of other solutions with finite 
contact angle for kq > kx. In this sense, kx is an eigenvalue for zero contact 
angle drops. 

From the numerical calculations we observe that the values of Sx are maxima 
for the solutions with zero contact angle. That means that for fixed V, m and 
7, the zero contact angle drops have maximum spreading speed (see Eq. I19|l . 

In Fig. 2 we show the solution with zero contact angle for several values of 
A = 1.1,1.5,2,3,4 and 5. The region where the slope tends to zero near the 
front is very small. There, the asymptotic expansion of the solution with zero 
contact angle is H{ri) = 

+lower order terms. 

In table 1 we give the values of Sx and Qx as a function of A. They may be 
useful for future experimental comparisons using Eq. 1191 and Eq. 1201 

Finally, within the lubrication approximation, the dissipation rate of energy 

is 

and for our self similar solution, it is proportional to 

2tt r i]\K'\^+^H^+'^dri. (23) 
Jo 

and then, using Eq. |^it is easy to see that the local rate of dissipation of energy 
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converges near the front. As a consequence, in contrast with a ncwtonian fluid, 
the paradox of the contact hne does nor arise for power law fluids with A > 1. 
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LIST OF FIGURES 

Fig 1- Solutions for A = 5/2 and k = 4, 3.00428 and 2 (from left to right). 

Fig 2- Zero contact angle solutions for A = 1.1, 1.5, 2.0, 3.0, 4.0 and 5.0 (from 
left to right). 
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